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ABSTRACT 

We have used the Sydney University Stellar Interferometer (SUSI) to measure the angular 
diameter of a Cir. This is the first detailed interferometric study of a rapidly oscillating A 
(roAp) star, a Cir being the brightest member of its class. We used the new and more accurate 
Hipparcos parallax to determine the radius to be 1.967 ± 0.066 Rq. We have constrained the 
bolometric flux from calibrated spectra to determine an effective temperature of 7420 ± 170 K. 
This is the first direct determination of the temperature of an roAp star. Our temperature is 
at the low end of previous estimates, which span over 1000 K and were based on either pho- 
tometric indices or spectroscopic methods. In addition, we have analysed two high-quality 
spectra of a Cir, obtained at different rotational phases and we find evidence for the presence 
of spots. In both spectra we find nearly solar abundances of C, O, Si, Ca and Fe, high abun- 
dance of Cr and Mn, while Co, Y, Nd and Eu are overabundant by about 1 dex. The results 
reported here provide important observational constraints for future studies of the atmospheric 
structure and pulsation of a Cir. 

Key words: stars: individual: a Cir - stars: fundamental parameters - stars: chemically pe- 
culiar - techniques: interferometric 
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1 INTRODUCTION 

The rapidly oscillating Ap (roAp) stars are chemically peculiar 
main-sequence stars with effective temperatures around 6500- 
8000 K. They are found around the classical instability strip but 
have quite different properties from the S Scuti stars, due to their 
slow rotation and strong global magnetic fields (typically of sev- 
eral kilogauss (kG), and up to 24 kG). They have peculiar abun- 
dances of Sr, Co, and certain rare earth elements such as Eu, Nd, 
Pr, Tb, an d Th that are found in h igh concentration high in the at- 
mosphere (Rvabchikova et al. 2007). The roAp stars present an in- 
triguing possibility to study element stratification, stellar evolution 
and pulsation in the presence of magnetic fields. 

The roAp stars oscillate in high overtone, low-degree p modes 
similar in period to the 5-minute acoustic oscillations in the Sun, 
but with coupling of the magnetic field, rotation and the oscilla- 
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tions that must be considered l lCunhal 20051) . The excitation of high- 
overtone pulsations, instead of the low overtones found in S Scuti 
stars, is thought to be directly related to the strong magnetic fields, 
which may suppress envelope convection, increasing the efficiency 
of the opacity mechanism in the r egion of hydrogen ionization 
jBalmforth et alJl200ll : ICunhall2002h . Asymptotic theory, valid for 
high overtones, predicts that regularly spaced peaks will dominate 
the frequency spectra of roAp stars, and this has b een observed in 
photo metric studies of several stars of the class ^Matthews et al.l 
1999). More recently, detailed pulsation studies of roAp stars have 
used high-resolution spectra collected at high cadence. From ana- 
lysis of the rare-earth element lines a complex picture has emerged, 
in which the mode amplitudes and sometimes phases depend on 
atmospheric hei ght, indicating the presence of running m agneto- 
acoustic waves l IKurtz et al . 2006; Rv abchikova e t al. 2007). 

a Circini (HR 5463, HD 128898; V = 3.2) is the brightest 
of 40 currently known roAp stars. It has one dominant oscillation 
mode with a period of 6.8 min and a semi-amplitude of 2.5 mmag 
in B. Four other modes with amplitudes lower by an order of mag- 
nitude w ere detected in an extensive three-site ground-based cam- 
paign bv lKurtz et al. Ul994l) . but a clear signature of the large sep- 
aration has so far not been found in the star. Recently, the WIRE 
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satellite observed a Cir which for the first time directly showed the 
rotational modulation (P ro t = 4.46 d) and a clear regular struc- 
ture of the pulsation peaks, with three modes having comparable 
amplitudes. The equidistant separatio n between them h as been in- 
terpreted as half the large separation (Bruntt et al. 2008). This new 
result offers the possibility to confront observations with theoreti- 
cal evolution and pulsation models for the star. 

To make progress in the asteroseismic modelling of a Cir we 
not only need to measure accurate frequencies with secure mode 
identification, but also to know with good accuracy the fundamen- 
tal parameters of the star. It is the goal of the current paper to 
determine the effective temperature, luminosity, mass and chem- 
ical composition of a Ci r. A core-wing ano maly is seen in the 
Balmer lines of Ap stars dCowlev et aifooil) . indicating that the 
temperature structure is different from normal A-type stars. There- 
fore, estimates of fundamental parameters from photometric in- 
dices o r spectral analysis are l ikely to be affected by systematic 
effects (Kochuk hov et alj|2005b . For example, T c g estimates for 
a Cir cover the range 7470 — 8730 K (considering la uncertain- 
ties; cf. Table |4](, spanning a large part of the range of T c s of the 
stars belonging to the roAp class. Also, esti mating the surface gr av- 
ity from a classical abundance analysis (e.g. Ku pka et al Jl 19960 by 
requiring lines of neut ral and ionized lines to y ield the same abun- 
dance is questionable (Rvabchi kova et ail2 002). 

The importance of obtaining accurate fundamental parameters 
through interfero metry for detailed ast eros eismic studies was dis- 
cussed in detail bv lCreevev et al.l (2007) and Cunha et al J 120071) . In 
the current paper we present new interferometric data to measure 
the angular diameter of a Cir. We use calibrated spectra to estimate 
the bolometric flux, which in turn allows us to determine the effec- 
tive temperature nearly independently of atmospheric models. In 
our compa nion paper on t he detection of the large frequency sep- 
aration jBruntt et al.ll2008l) . we apply the parameters found here in 
our theoretical modelling of the star. 

The presence of strong magnetic fields in roAp stars is thought 
to su ppress the turbulence in their outer atmospheres dMichaud 
1970). For this reason, enhanced d iffusion leads to stratification of 
elements (e.g. iBabel & Lanzlll992l) . A more complicated pictur e 
was found in the roAp star HR 3831 bv lKochukhov et al.l J2004) . 
who monitored significant changes in spectral lines with rotational 
phase. Using a Doppler imaging technique they found rings and 
spots of different spatial scales, with enhanced ab undances up to 7 
dex for some elements (see also Liift inger et al.l2007l) . Such studies 
require high-quality data covering a wide range of rotational phases 
and this has not been done for a Cir. Single spectra indicate that it is 
a typical roAp star, with high abu ndances of certain elements such 
as Co and Nd dKupka et al.ll996|). To confirm the evid ence of spots 
on q Cir found by Kochuk hov & Rvabchikoval bOOll) . we present 
in Sect.[6]an analysis of two high-resolution and high-S/N spectra 
taken at different rotational phases. 



2 OBSERVATIONS AND DATA REDUCTION 



Table 1. Adopted parameters of the two calibrator stars. The primary cali- 
brator, f3 Cir, was used on all nights while the secondary calibrator, a Lup, 
was also used on 22 April and 30 July. 



The S ydney University Stellar Interferometer (SUSI; iDavis et al] 
1999) was used to measure the squared visibility, i.e. the nor- 
malised squared modulus of the complex visibility or V 2 , on a 
total of eight nights. The red- table beam-combination system was 
employed with a filter of centre wavelength and full-width half- 
maximum 700 nm and 80 nm, respectively. This system, including 
the standard SUSI obs erving procedure , data reduction and calibra- 
tion, is described bv lDavisetaT] ( l2007l) . 



HR 


Name 


Spectral 


V 


UD Diameter 


Separation 






Type 


(mag) 


(mas) 


from a Cir 


5670 


/? Cu- 


A3 V 


4.05 


0.58 ± 0.05 


7? 42 


5425 


ff Lup 


B2III 


4.42 


0.20 ± 0.04 


14? 59 



Target observations were bracketed with calibration measure- 
ments of nearby stars, which were chosen to be essentially un- 
resolved. Using an intrinsic colour interpolation (and spread in 
data) of measur ements made with the Na rrabri Stellar Intensity 
Interferometer dHanburv Brown et al.ll 19740 , the angular diameters 
(and associated uncertainty) of the calibrator stars were estimated 
and corrected for the effects of limb-darkening. However, post- 
processing indicated that two of the calibrator stars (HR 4773: 
7 Mus and HR 5132: e Cen) are binaries and could not provide ad- 
equate calibration information. The adopted stellar parameters of 
the remaining calibrator stars are given in Table[T] A linear interpo- 
lation of the calibrator transfer functions was then used to scale the 
observed squared-visibility of the target. This procedure resulted 
in a total of 60 estimations of V 2 and a summary of each night is 
given in Table [2] 

This calibration is not ideal due to (i) the latency between cal- 
ibration and target measurements; and (ii) the use of only one cali- 
brator on all but two nights. However, the use of a nonlinear fu nc- 
tion to correct partially for residual seeing effects ( Ireland 2006) as 
part of the data reduction reduces any systematic error introduced. 

It should also be noted that the uncertainty in the calibrator 
angular diameters only has a small effect on the final, calibrated 
target V 2 . Since our calibrators are nearly unresolved, their true 
visibilities are close to unity and depend only weakly on their actual 
diameters. Therefore, propagation of these uncertainties is greatly 
reduced and only has a minor effect on the final uncertainty of the 
diameter of a Cir. 



3 ANGULAR DIAMETER 

In the simplest approximation, the brightness distribution of a star 
can be modelled as a disc of uniform (UD) irradiance with angu- 
lar diameter #ud- A two-aperture interferometer has a theoretical 
response to such a model given by 



\v\ 2 = 



A 



2JiO|6|0/A) 



7r|fe|6»/A 



(1) 



where A is the observing wavelength, b is the baseline vector pro- 
jected onto the plane of the sky and Ji is a first-order Bessel 
function. For stars with a compact atmosphere, only small correc- 
tions are needed t o account for monochromatic limb-darkening (see 
IDavis et al.l2 000). The parameter A is included in the model to ac- 
count for (i) any incoherent flux from the nearby faint K5 V star 
GJ 560 B that could enter the interferometer's field of vie and 
(ii) any instrumental effects arising from the differing spectral types 
of a Cir and the calibrator stars. The latter effect was found to be 



iPoveda et al ] ll994l) list the magnitude V = 8.47, while 
iGould & Chaname 1 2004) list V = 9.46 and a separation from a Cir of 
15.7". 
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Table 2. Summary of observational interferometry data for a Cir from 
SUSI. The night of the observation is given in Columns 1 and 2 as a cal- 
endar date and the mean MJD— 54 000. The nominal and mean projected 
baselines in units of metres are given in Columns 3 and 4, respectively. 
The weighted-mean squared visibility, associated error and the number of 
observations during each night are given in the last three columns. 



Date 




Nominal 


Projected 








in 2007 


MJD 


Baseline 


Baseline 


V 2 




N 


April 18 


208.62 


5 


4.08 


0.964 


0.012 


10 


April 20 


210.63 


20 


16.37 


1.039 


0.016 


9 


April 22 


212.66 


20 


16.29 


0.900 


0.015 


9 


May 24 


244.56 


80 


65.25 


0.533 


0.008 


8 


May 25 


245.55 


80 


65.34 


0.528 


0.010 


8 


June 22 


273.46 


5 


4.10 


0.987 


0.024 


4 


July 30 


311.40 


10 


32.48 


0.845 


0.038 


6 


August 4 


316.40 


80 


64.42 


0.560 


0.013 


6 



negligible as the spectral types are very similar (see also lDavis et al.l 



neglig 



Eq. [T] is strictly valid only for monochromatic observations, 
but is an excellent appr oximation when the e f fective wavelengt h 
is correctly defined (see Tango & Davis 20021: iDavis et al]|2007l) . 
During the a Cir observations, the coherent field-of-view of SUSI 
was found to be greater than 6.6 mas, which is larger than the an- 
gular extent of a Cir (see below). Hence, bandwidth smearing can 
be considered negligible. The interferometer's effective wavelength 
when observing an A-type main -sequence star is approximately 
696.0 ± 2.0 nm jDavis et al .1120071) . 

We estimated (9ud and A by fitting Eq. [T] to the 60 individ- 
ual measures of V 2 . The fit was achieved using a \ 2 minimisation 
with an implementation of the Levenberg-Marquardt method. The 
reduced \ 2 of the fit was 2.2, implying that the measurement un- 
certainties were underestimated. We have therefore multiplied the 
V 2 measurement uncertainties by f x = \/2.2 to obtain a reduced 
X 2 of unity. In Fig. Q] the fitted uniform disc model (solid line) is 
compared to the weighted-mean V 2 measures at each of the four 
baselines. The bottom panels show the individual measures at each 
baseline. The plotted error bars are scaled by f x . 

Formal uncertainties in 9\jd and A, derived from the diago- 
nal elements of the covariance matrix, may be underestimates for 
three reasons: the visibility measurement errors may not strictly 
conform to a normal distribution, Eq. Q]is non-linear, and the pri- 
mary calibrator is partially resolved. The effect of these factors on 
the model parameter uncertainties was investigated using Markov 
chain Monte Carlo (MCMC) simulations. This method involves a 
likelihood-based random walk through parameter space, using an 
implementation of a Metropolis-Hastings algorithm, and yields the 
full marginal posterior probability densit y function (PD F; for an 
introduction to MCMC see Chapter 12 of lGregory||2005h . Further- 
more, the current knowledge of the system can be included in the 
analysis by assuming an a priori distribution. For example, knowl- 
edge of the effective wavelength and calibrator angular diameter 
(and associated uncertainties) can be included into the uncertainty 
estimation of the model parameters. 

Over 25 MCMC simulations were completed, each with 10 6 
iterations, and adopting Gaussian likelihood distributions for the ef- 
fective wavelength and calibrator angular diameters. The resulting 
PDFs for #ud produced a la parameter uncertainty approximately 
twice as large as the formal uncertainty derived from the covari- 
ance matrix, while those for A were approximately equal. Since 




16.0 16.4 32.0 32.8 63.5 
Projected baseline (m) 

Figure 1. The top panel shows the baseline-averaged V 2 measures (open 
circles) with the fitted uniform disc model overlaid (solid line). The bottom 
panels show the individual V 2 measures (filled circles) at each baseline. 



the simulations take into account the uncertainties in the effective 
wavelength and calibrator angular diameter, we adopt the standard 
deviations from the MCMC simulations. We believe they represent 
the most realistic and conservative parameter uncertainty estimates 
for our dataset. 

Final values oftfuD and A are 1.063±0.034 masand0.992± 
0.006, respectively. The value of A, while close to unity, indicates 
that the measurements are affected slightly by the incoherent flux 
of the faint star GJ 560 B. 

The work of IDavis et al] J2000l) can be used to correct the 
angular diameter for the effect of limb-darkening, provided the 
effective temperature (T e g), surface gravity (log g) and metallic- 
ity ([Fe/H]) are known. Since roAp stars have non-standard atmo- 
spheres and lDavis et al. (2000) used ATLAS9 atmospheric models, 
we investigated the limb-darkening correction within the following 
large region of parameter space: T e g = 7250 to 8500 K, log g — 
4. to 4.5, [Fe / H] = —0.3 to +0.3. The limb-darkening correction 
of IDavis et a at an observing wavelength 696 nm, varies 

between 1.0344 and 1.0426 and is more sensitive to changes in 
T e g than in log<? or [Fe/H]. We conservatively set an uncertainty 
of 0.010 in the limb-darkening correction, but note that this un- 
certainty is small compared to the uniform-disc angular diameter 
uncertainty. 

For the fundamental parameters found in Sect. [5] we get a 
limb-darkening correction of 1.039 ± 0.010, to obtain the limb- 
darkened (LD) angular diameter of 6>ld = 1.105 ± 0.037 mas. 



4 THE BOLOMETRIC FLUX OF a Cir 

We derived the bolometric flux of a Cir using two different 
datasets. The first value was obtained by combining the observed 
ultraviolet flux retrieved from IUE Newly Extracted Spectra (INES) 
data archive, with the theoretical flux obtained from a Kurucz 
model (with IDL routine KURGET1). T o compute the integrated 
fluxes, we followed the same method as INorthI dl98lh . The IUE 
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Table 3. Physical parameters of a C ir. All parameters a re from the current 
study except the Hipparcos parallax Ivan Leeuwen 2007). 



Parameter 


Value 


Uncertainty (%) 


6» LD (mas) 


1.105±0.037 


3.4 


/bol(10- 9 Wm- 2 ) 


1.23±0.07 


5.7 


7T P (mas) 


60.36±0.14 


0.2 


T cB (K) 


7420±170 


2.3 


R (R©) 


1.967±0.066 


3.4 


L(Lq) 


10.51±0.60 


5.7 


M(M Q ) 


1.7±0.2 


12 



measurements were used to compute the flux in the wavelength in- 
terval 1150 A < A < 3207 A, while the flux from A = 3210 A to 
infinity was computed from the Kurucz model that best fitted the 
seven intrinsic colours of the star i n the Geneva p hotometric sys- 
tem, retrieved from the catalogue o f lRufenerl fl989). Moreover, the 
energy distribution was extrapolated to the interval 912 A < A < 
1150 A, assuming zero flux at A = 912 A. The resulting integrated 
flux was /boi = (1.28 ± 0.02) x 10~ 6 erg cm " 2 s"\ where the 
quoted uncertainty range corresponds to the difference between the 
maximum and minimum integrated fluxes allowed by the uncer- 
tainties in the parameters of the Kurucz model adopted, which, in 
turn, reflect the uncertainties in the photometric data. 

The abnormal flux distributions that are characteristic of Ap 
stars make determinations such as that described above, based 
on atmospheric models appropriate to normal stars, rather unreli- 
able. When spectra calibrated in flux are also available for visi- 
ble wavelengths, a more reliable determination of the bolometric 
flux can in principle be obtained. Unfortunately, even though two 
low-resolution spectra calibrated in flux are available in the litera- 
ture f or a Cir (see catalogues by Alekseeva et al1ll996l : lBurnashevl 
1 1985b . the errors associated with the calibrations in flux as func- 
tion of wavelength are not given in the corresponding bibliographic 
sources. Nevertheless, we have calculated a second value for the 
bolometric flux using the same method as above, but replacing the 
synthetic spectra obtained for the Kurucz model by the two low- 
resolution spectra of a Cir, calibrated in flux, retrieved from the cat- 
alogues mentioned above. In this case, the flux as function of wave- 
length was obtained by combining the IUE measurements in the 
wavelength interval 1150 A < A < 3349 A, the mean of the two 
low-resolution spectra calibrated in flux in the interval 3200 A < 
A < 7350 A, and the flux derived from the mean of the two Kurucz 
models that best fitted each of the two low-resolution spectra, for 
wavelengths longer than 7370 A. The integrated flux obtained in 
this way was /boi = (1.18±0.01) X 10 _6 erg cm ~ 2 s _1 , where the 
uncertainties reflect the difference between the integrated fluxes de- 
rived when each of the two low-resolution spectra (and correspond- 
ing best-fit Kurucz model at higher wavelengths) were considered 
separately. 

Since the details and associated uncertainties of the calibra- 
tions in flux of the low-resolution spectra used in the latter cal- 
culation are unknown, we will take a conservative approach and 
allow the bolometric flux to vary within the two extremes of the 
values derived through the two different approaches. Hence, in 
what follows, we use for the bolometric flux, the mean value 
/boi = (1.23 ± 0.07) x 10 _6 erg cm " 2 s -1 . 
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Figure 2. The position of a Cir in the Hertzsprung-Russell diagram, with 
three evolution tracks for models with masses of 1.70, 1.715 and 1.73 Mq 
for comparison. The constraints on the fundamental parameters are indi- 
cated by the lcr-error box (T e ff, L/Lq) and the diagonal lines (radius). 



5 FUNDAMENTAL PARAMETERS OF a Cir 

We use the measured angular diameter and bolometric flux to esti- 
mate T e g from the definition 



CrT e 4 ff = 4/ bo l / 6» L 



(2) 



where a is the Stefan-Boltzmann constant. We can also obtain es- 
timates of the radius and luminosity from the measured angular 
diameter and bolometric flux, when they are combined with the 
parallax: 



R = 
and 



26C/n p , 



(3) 



(4) 



L = 4-7r/ bo i C /tt p , 

where C is the conversion from parsecs to metres. We have used 
a parallax of n p — 60.36 ± 0. 14 mas, based on th e new reduc- 
tion of the raw Hipparcos data ( Ivan Le euwen 2007). It has lower 
uncertainty, but is in good agreement with the original data release 
dES All 19971) . which was tt p = 60.97 ± 0.58 mas. The values we 
have determined for the fundamental parameters of a Cir, along 
with the observable quantities, are given in Table [3] Due to the ac- 
curate parallax, the uncertainties in the radius and luminosity are 
dominated by the uncertainty in the angular diameter and bolomet- 
ric flux, respectively. The uncertainty in the effective temperature 
is determined by the uncertainties in both the angular diameter and 
the bolometric flux. 

In Fig. [2] we show the location of a Cir in the Hertzsprung- 
Russell diagram. The la-error box and diagonal lines are based 
on the fundamental parameters and asso ciated errors derived in 
this study (see Table 0. Three CESAM dMorellll997h evolution- 
ary tracks crossing the error box for a Cir are also shown. The 
tracks were calculated for solar metallicity using a mixing length 
parameter of a = 1.6. Diffusion and rotation were not included. 
Note that changes in the model input parameters can significantly 
modify the evolutionary tracks for models with the sa me mass. In 
a study of the roAp star HR 1217. ICunha et al.l j2fJ03h found that, 
while changes in the mixing length parameter and convective over- 
shoot had a minor impact on the mass derived from model fitting 
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Table 4. Fundamental atmospheric parameters of a Cir as found in the literature and the current study. 



Study 


T eff [K] 


log 3 


[Fe/H] 


Applied method 


Kurtz & Martinez (1993) 


8000 






Hp index 


North et al. (1994) 


8000 


4.42 


+0.2 


Geneva indices 


Kunka et al. (1996) 


7900±200 


4.2 ±0.15 


-0.13 ±0.10 


Spectral analysis 




7950 


4.3 




Stromgren indices 




7880 


4.4 


+0.3 


Geneva indices 


Sokolov (1998) 


8440±290 






Continuum slope 




7680 






Geneva indices 


Matthews et al. (1999) 


8000±100 






Hp index 


Kochukhov & Bagnulo (2006) 


7673±200 


4.12±0.06 




Geneva indices 


This study 


7420±170 


4.09±0.08 




Interferometry + bolometric flux 








-0.12 + 0.08 


Spectral analysis 



of the classical observables, the uncertainties in the initial helium 
abundance and in the global metallicity of the star resulted in a sig- 
nificant uncertainty in the mass determination that in the worst case 
could be as large a s 10%. Therefore, from Fig. [2] and the results of 
ICunha et all (2003). we conservatively estimate the mass of a Cir 
to be M = 1.7 ± 0.2 Mq. Even so, we obtain an accurate sur- 
face gravity of logg = 4.09 ± 0.08. This is a consequence of the 
accurately measured interferometric radius. 



5.1 Comparison with previous estimates 

Previous studies of a Cir have used either spectroscopy or pho- 
tometric indices to estimate the fundamental atmospheric param- 
eters. _We_list_sj5me_receri^stim in Ta- 
blel4l lKurtz & Martinez! d 19931) and lMatthews et all d 19991) used the 
Stromgren Hp index to estimate r P ff. lNorth et alj d 19941) made an 
abundance analysis of a Cir based on a few selected wavelength 
regions with high resolution, but the fundamental param eters were 
fixed based on Geneva photometry. iKupka et all d 19961) analysed 
high-resolution spectra of a Cir and estimated the fundamental 
parameters using Geneva and Stromgren photometric indices, the 
Balm er lines and a de tailed classical analysis of several iron-peak 
lines. ISokolovl {1998) proposed a new method to estimate T e g of 
roAp stars based on the slope of the conti nuum near the Balmer 
jump and found a high T e g for a Cir. iKochukhov & Bagnulol 
(2006) analysed a large sample of magnetic chemically peculiar 
stars and they used Geneva indices to estimate the fundamental 
parameters of a Cir. Some of the studies mentioned here that 
used photometric indices have not given uncertainties on the pa- 
rameters. Realistic values are about 200 K, 0.2 dex, and 0.1 de x 
on T eff , logfl, and [F e/H], respectively dKupka & Brunttl 1200 ll) . 
IKochukhov et all d2005l) explored the influence of magnetic fields 
on the photometric indices and found relatively small corrections, 
especially for field strengths of a few kG or less, as found in a Cir 
(see Sect.|6~TV 

Previous estimates of \ogg lie in the range from 4.1 to 4.6, 
considering la error bars. We have obtained an accurate estimate, 
log g= 4.09 ± 0.08, which is lower than previous estimates but 
in general agreement. The range in T c g estimated for a Cir in the 
literature is 7470 to 8730 K, considering la error bars. This has 
hampered detailed studies of a Cir that rely on theoretical mod- 
els, such as spectroscopic analysis (atmospheric models) or as- 
teroseismology (evolution and pulsation models). In general the 
Geneva colours give lower T c h than the Stromgren or Hp in- 
dices. The slight differences in the parameters found using the same 
photometric system are due to different calibrations being used. 



Table 5. The atomic number, element name, wavelength, and oscillator 
strength (log gf ) from the VALD database for lines used in the abundance 
analysis. The full version of this Table is available in the on-line version. 



El. 


A [A] 


log af 


El. 


A [A] 


logg/ 


El. 


A [A] 


log gf 


6 Ci 


5017.090 


—2 


500 


Till 


5490.690 


-2 


650 


Fei 


5217.389 


-1.070 


Ci 


5540.751 


— 2 


376 


23 VI 


5670.853 


-0 


420 


Fei 


5232.940 


-0.058 


Ci 


5603.724 


—2 


418 


Vi 


6081.441 


-0 


579 


Fei 


5242.491 


-0.967 


Ci 


5643.368 


— 2 


670 


VI 


6199.197 


-1 


300 


Fei 


5250.646 


-2.181 


Ci 


5800.602 


-2 


338 


24 Crl 


5204.506 


-0 


208 


Fei 


5253.462 


-1.573 


Ci 


6413.547 


— 2 


001 


Cn 


5224.972 


— 


096 


Fei 


5263.306 


-0.879 


»Oi 


5330.737 


-0 


984 


Cn 


5247.566 


-1 


640 


Fei 


5269.537 


-1.321 


Oi 


6155.986 


— 1 


120 


Cn 


5296.691 


— 1 


400 


Fei 


5281.790 


-0.834 


Oi 


6156.776 


-0 


694 


Cn 


5297.376 


+ 0.167 


Fei 


5283.621 


-0.432 


Oi 


6158.186 


-0 


409 


Cn 


5348.312 


— 1 


290 


Fei 


5288.525 


-1.508 


11 Na I 


5682.633 


-0 


700 


Cn 


5787.965 


-0 


083 


Fei 


5302.302 


-0.720 


Nal 


5688.205 


-0 


450 


Cn 


6661.078 


-0 


190 


Fei 


5315.070 


-1.550 


12 Mg. 


5172.684 


-0 


402 


Cm 


5237.329 


-1 


160 


Fei 


5324.179 


-0.103 


Mgl 


5183.604 


-0 


180 


Cru 


5246.768 


— 2 


466 


Fei 


5353.374 


-0.840 


Mgl 


5528.405 


-0 


620 


Cm 


5305.853 


-2 


357 


Fei 


5364.871 


+ 0.228 


Mgl 


5711.088 


— 1 


833 


Cm 


5308.408 


— 1 


846 


Fei 


5365.399 


-1.020 


^All 


5557.063 


-2 


110 


Cm 


5310.687 


-2 


280 


Fei 


5367.467 


+ 0.443 


All 


6696.023 


— 1 


347 


Cm 


5313.563 


— 1 


650 


Fei 


5373.709 


-0.860 


14 Sii 


5665.555 


-2 


040 


Cm 


5334.869 


-1 


562 


Fei 


5383.369 


+ 0.645 


Si I 


5708.400 


— 1 


470 


Cm 


5407.604 


— 2 


151 


Fei 


5386.334 


-1.770 


Si I 


5772.146 


-1 


750 


Cm 


5420.922 


-2 


458 


Fei 


5389.479 


-0.410 


Si I 


6125.021 


-0 


930 


Cm 


5478.365 


— 1 


908 


Fei 


5391.461 


-0.825 


Si I 


6142.483 


-0 


920 


Cm 


5508.606 


-2 


110 


Fei 


5393.168 


-0.715 


si i 


6414.980 


— 1 


100 


Cm 


5510.702 


— 2 


452 


Fei 


5398.279 


-0.670 


Si I 


6856.037 


— 1 


060 


Cm 


6068.023 


-1 


736 


Fei 


5400.502 


-0.160 


16 Si 


6395.163 


— 2 


040 


Cm 


6070.105 


-2 


944 


Fei 


5405.775 


-1.844 


Si 


6396.644 


— 1 


820 


25 Mni 


5377.637 


-0 


109 


Fei 


5410.910 


+ 0.398 


Si 


6415.522 


— 1 


360 


Mm 


5413.689 


-0 


589 


Fei 


5415.199 


+ 0.642 


Si 


6743.531 


-0 


9 20 


M ii I 


5420.355 


-1 


462 


Fei 


5424.068 


+ 0.520 


Si 


6757.171 


-0 


310 


Mn I 


5470.637 


-1 


702 


Fei 


5434.524 


— 2.122 


20 Ou 


5512.980 


— 


712 


M ii I 


5537.760 


-2 


017 


Fei 


5466.396 


— 0.630 



Our value, T cff = 7420 ± 170 K, is the lowest of all temperature 
determinations. However, it is in acceptable agreement with the 
T e ff from Geneva indices (the mean value from four estimates is 
78 10 + 160 K) and ba rely agrees with T cff = 7900 ± 200 K found 
by iKupkaetaild 19961) based on spectroscopic analysis of Fe I and 
Fe II lines. We note that a higher T ff is also implied in our anal- 
ysis (in Sect. [6j of Fe I spectral lines due to a strong correlation 
with excitation potential. This could either imply stratification of 
Fe, as seen in other roAp stars, or, alternatively, that T c s of the at- 
mospher ic model is too low by 600 K. We believe the higher T e a 
found by Kupka et al. ( 1996) is explained by this. 

We stress that our estimates of logg and T c g are largely 
model-independent, unlike the previous estimates. Therefore we 
recommend that these values be used in future studies of a Cir. Our 
value for T a g depends on the angular diameter from interferometry 
and the bolometric flux. The latter estimate should be confirmed by 
making new observations. 
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6 SPECTRAL ANALYSIS OF a Cir 

We have obtained high quality spectra of a Cir with the Ultraviolet 
and Visual Echelle Spectrograph (UVES) on the Very Large Tele- 
scope (VLT). The s pectra were colle cted as part of asteroseismic 
campaigns in 2004 dKurtz et alj |2006) and 2006. In the latter case, 
the spectra were used to compare the variability of a Cir in radial 
velocity and photometry f rom the star tracke r on the WIRE satel- 
lite and from the ground teruntt et al.ll2008l) . The spectrum from 
2004 is comprised of 64 co-added spectra collected over 0.5 h, at 
a mean rotational phase of 0.230. We define phase zero to be at 
minim um light in the r otational modulation curve measured with 
WIRE iBruntt et al]2008h . The spectrum from 2006 is comprised of 
149 co-added spectra collected over 2.4 h at mean rotational phase 
0.975, i.e. close to minimum light. The resulting spectra have S/N 
of 800 and 1100, respectively, in the continuum around 5500 A. 
The spectral resolution is about R = 110 000 and we used spec- 
tra from the red arm in the wavelength ranges 4970 — 5950 A and 
6090 - 6900 A. We avoided the range 5880 - 5950 A, which is 
affected by many telluric lines. This was especially the case for the 
2006 spectrum, which was collected through cloud. For this reason, 
we base our final results on the spectrum from 2004. The spectra 
were carefully normalized by identifying continuum windows in a 
synthetic spectrum with the same parameters as a Cir. 

The spectral line s were analysed using the VWA pack- 
age l lBruntt et al.ll2004l l2008h . The abundance analysis relies on 
the calculation of s ynthet ic spectra with the SYNTH code by 
IValenti & Piskunovl (1996). Atomic parameters and line broad- 
ening coefficients were extracted from the VALD database 
jKupkaetal.il 19991) . Table [5] lists the lines we have used in the 
analys is. We used modified ATLAS9 models from iHeiter etldl 
(2002) with the adopted parameters from interferometry T c g = 
7420 ± 170 K and log g = 4.09 ± 0.08. We initially assumed solar 
metallicity, but adjusted the abundances of individual elements in 
the model after a few iterations. We adjusted the microturbulence 
to £ t = 1.60 ± 0.15kms _1 by requiring that there is no correla- 
tion between the abundance and strength of 50 weak Fe I lines with 
equivalent widths < 90 mA and excitation potentials in the range 
3-5 e V. We repeated our a nalysis with VWA using the LLModels 
code JShulvak et afl l2004). and found the same result for all ele- 
ments, with differences below 0.03 dex. 

The abundances relative to the Sun for 23 specie jf] are given 
in Table [6] To estimate uncertainties in the abundances we mea- 
sured abundances for model atmospheres with higher values of T e g 
(+300 K), logg (+0.3 dex), and microturbulence (+0.4kms -1 ), 
and made a linear interpolation using the uncertainties of these pa- 
rameters. In addition, we added Q.YJ /\fN dex to the uncertainty, 
where N is the number of lines, and 0.17 is the average rms stan- 
dard deviation on the abundances for four elements with the most 
spectral lines, i.e. Ca, Fe, Co, and Nd. For example, for Fel the 
estimated uncertainty on the abundance is 0.08 dex, including con- 
tributions from the uncertainty on the fundamental parameters. 

The Fel lines are the most numerous in the observed spec- 
tra. We find a significantly lower abundance from Fe I lines with 
low excitation potential, e. For lines with 1 < e [eV] < 3 the 
abundance relative to the Sun is —0.33 ± 0.02 (28 lines) and for 
3 < e [eV] < 5 the abundance is -0.07 ± 0.02 (89 lines). The 
quoted uncertainties are the rms on the mean value, and represent 

2 We quote the abundance of an ion A as [A] = logA'A/A'tot — 
(log N\ /Nto t ) p, , where N i s the number of atoms and the solar values 
are taken from Grevesse e t alj |2007|) . 



Table 6. Abundance of a Cir relative to the Sun based on the analysis of the 
spectrum from 2004. The mean abundance and the number of neutral and 
ionized lines are denoted I and II. 



Element 


[A]i 


Ni 


[A] n 


Nil 


C 


+0.05 ± 0.07 


6 






o 


-0.15 ± 0.09 


4 






Na 


-0.21 ± 0.14 


2 






Mg 


-0.01 ± 0.17 


2 






Al 


+0.18 ± 0.19 


2 






Si 


-0.14 + 0.20 


7 






s 


+0.35 ± 0.08 


5 






Ca 


+0.19 ± 0.10 


20 






Sc 






-0.42 + 0.11 


3 


Ti 


-0.20 + 0.13 


5 


-0.17 + 0.09 


5 


V 


+0.60 + 0.21 


3 






Cr 


+0.48 ± 0.11 


8 


+0.49 + 0.10 


11 


Mn 


+0.68 ± 0.14 


7 


+0.68 + 0.09 


4 


Fe 


-0.12 ± 0.08 


125 


-0.03 + 0.05 


20 


Co 


+1.20 + 0.08 


21 






Ni 


-0.58 ± 0.10 


7 






Cu 


-1.02 ± 0.18 


2 






Y 






+0.88 + 0.12 


10 


Zr 






+0.37 + 0.23 


1 


Ba 






-0.70 + 0.18 


2 


Ce 






+0.74 + 0.13 


3 


Nd* 






+1.03 + 0.10 


12 


En 






+1.58 + 0.15 


2 



*From three Nd III lines we find [Nd]= +2.06 ± 0.12. 



only the internal error. The difference between low and high ex- 
citation lines is significant and could be an indication of vertical 
stratification, meaning the abundance of Fe is higher in the deeper 
layers, as discussed by Rva bchikova etalj ( l200ll2007h . 

In Fig.[3]we show a small region of the observed spectra from 
2004 (solid thick line) and 2006 (dashed thick line), with identi- 
fication of the elements causing the lines. We show two synthetic 
spectra (thin lines) with different abundances of cobalt and the rare- 
earth element neodymium. In the first synthetic spectrum the abun- 
dances are solar and in the second we have adjusted the abundances 
of Co and Nd to fit the observed spectrum from 2004. We find over- 
abundances of Col and Ndll of +1.20 ± 0.08 and +1.03 ± 0.11 
using 21 and 12 line s, respectively. This is in acceptable agreement 
with the analysis bv lKupka et id] jl996h who found +1.62 + 0.20 
and +1.24 + 0.25, using 6 lines of each ion. For other ions listed in 
Table[6]we find genera l agreement with the previous detailed study 
bv lKupkaetai] ( ll996h . 

The two spectra we analysed yield the same mean abundances, 
when averaging over several lines, with the mean abundances dif- 
fering by less than 0.05-0.10 dex in the two spectra. A difference 
of 0.20 dex was found for Vl and EuII, but the analysis is based 
on just 3 and 2 lines, respectively. It is interesting to note that the 
equivalent widths and detailed shapes of the Co I and Nd II lines are 
significantly different in the two observed UVES spectra shown in 
Fig. [3] This indicates the presence of sp o ts on the surface, as first 
reported by Kochukho v & Rvabchikov3 ( 1200 ll) . A more detailed 
study of a Cir is required to explore this, as was done for the roAp 
star HR 3831 bv lKochukhov et all J2004I) . 
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Figure 3. The two observed UVES spectra of a Cir are compared (solid and dashed thick lines) in a region with strong neodymium and cobalt lines. Also 
shown are two synthetic spectra (thin lines) for solar abundance and when the abundance of Nd and Co is increased by about a factor of 10. 



6.1 Limits on the magnetic field strength 

The roAp stars are magne tic, with surface fi eld strengths as high 
as 24.5 kG in HD 154708 faubrig et alj|2005h . but more typically 
of several kG fsee lKurtz et al.ll2006T) . Knowledge of the magnetic 
field strength and structure is important for the study of a Cir. 
A variable longitudinal magn etic fiel d was discovered in t his sta r 
by iBorra & Landstreetl d 19751) and IWood & Campusand Jl975h . 
The latter paper shows a variable field with negative polarity and 
field strength extremum of — 1.5kG, but with a large uncertainty. 
Borra & Landstreetl dl 9751) obtained more precise results using a 
photoelectric B aim er line polarimeter and could not confirm the ex- 
tremum obtained bv lWood & Campusand dl975t) . They made seven 
measurements over seven consecutive nights, from which they de- 
tected a variable field ranging from —560 to +420 G. From these 
observations they proposed a rotation period of the o rder of 12 d - 
a value we now know to be incorrect. iMathyj < fl99lt) attempted to 
test the possibility of a short rotation period close to 1 d by observ- 
ing the longitudinal field strength. He obtained five spectra in three 
nights and di d not detect any magnetic field within the errors of his 
observations. iMathvs & Hubrid dl997h presented two more longi- 
tudinal field measurements which show a negative polarity mag- 
netic field, but their results were still below the 3a level. 

From these observations it is clear that the star does have 
a magnetic field and that the longitudinal component is rather 
weak. The precision of the published observations is insuffi- 
cient, however, to prove variability with the known rotation pe - 
riod of P rot = 4 46 d dKurtz et al.1 1 1994 iBruntt et al.1 2008). 
Bvch kov et alj d2005l) . for example, combined the data and derived 
a magnetic variation period equal to the known rotation period, but 
we conservatively are not convinced that rotational magnetic vari- 
ability has been found, given the error bars of the individual mea- 
surements. While there is an indication of a rotational magnetic 
curve, additional high-precision observations of the longitudinal 
magnetic field for this bright star are sorely needed. 

Estimating the magnetic field modulus (H) for a Cir is also 
not a simple task, given its weak longitudinal field and rather high 
(for an roAp star) rotational velocity of wsini = 13kms _1 . No 
Zeeman components can be seen directly in our high resolution 
spectra, so we searched for the effect of magnetic broadening by 
measuring the full-width at half maximum (FWHM) for selected 
unblended spectral lines with a spread of Lande factors. This gives 
a value of the magnetic field similar to the mean magnetic field 



modulus obtained from split Zeeman components in slowly rotat- 
ing stars where the components are resolved. From the measure- 
ment of the FWHM for 20 Fe and Cr lines we obtained a magnetic 
field modulus of 0.9 ± 0.5 k G, a significant imp rovement on the 
upper limit of 3-3.5 kG set bv lKupka et all (l996), using a similar 
approach. 

We estimated an upper limit for (H) by comparing average 
observed spectral lines with calculated lines in the presence of a 
magnetic field using the SYNTHMAG code jPiskunov|[l999l) . The 
line of Cm 5116 A with a large Lande factor of g = 2.92 is rather 
weak and asymmetric with a stronger blue part of the profile in 
the first observing set of 2004 and with a stronger red part for the 
second observing set from 2006. The synthetic line profile fit to this 
line using SYNTHMAG needed no magnetic field strength greater 
than 1-1.5 kG. Line profiles with large Lande factors belonging 
to iron, such as Fe I at 5501 A, 5507 A and 6337 A, fitted slightly 
better for a magnetic field strength up to 2 kG, in comparison with 
synthetic spectra with zero magnetic field. We therefore place an 
upp er limit of (H) < 2 kG. 

iRomanvukl |2004) compared longitudinal magnetic field ex- 
trema with the mean magnetic field modulus for 39 Ap stars and 
found a linear relation: (H) — 1.01 + 3.16(Z/;} kG. Using this re- 
lation an d the extremum of the longi tudin al field measurements ob- 
tained bv lBorra & Landstreetl d 1975b and lMathvs & Hubrid d 19971) 
gives a mean magnetic field modulus of 1.8 kG and 1.3 kG, respec- 
tively. These values are consistent with the upper limit estimated 
from our UVES spectra. 



6.2 Simplifications: magnetic fields and NLTE 

All atmospheric models include simplifications to the physical con- 
ditions in real atmospheres. This is especially the case for the roAp 
stars, in which magnetic fie lds may alter the atmospheric struc- 
ture dKochukhov et alj 120051) and where obse rvational evidence 
for st ratification of certain elements is found dRvabchiko va et ail 
2002). More realistic models of magnetic stars are being devel- 
oped dKochukhov et alj|2005l) , but they indicate that the structure 
in a Cir should only be mildly affected by the relatively weak mag- 
netic field of ^ 2 kG that we determined above. 

The ATLAS9 models we applied assume local thermody- 
namical equilibrium (LTE), but at T c g above 7000 K non-LTE 
effects become important. For a normal star with solar metal- 
licity and T e ff = 7500 K, the correction for neutral iron is 
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[Fe I/H]nlte = [Fe I/H]lte+0.1 dex jRentzsch- Holml 19961) . This 
seems to agree with the slight difference in Fe abundance we find 
from neutral and ionized lines. We adopt this correction for our fi- 
nal result for the photospheric metallicity, which is the mean of the 
125 Fe I lines in a Cir and, hence, [Fe/H] = -0.12 ±0.08. We note 
that the overall metallicity of a Cir, below its photosphere, is much 
more uncertain. 



7 CONCLUSION 

We present the first angular diameter measurement of an roAp star, 
a Cir, based on interferometric measurements with SUSI. In com- 
bination with our estimate for the bolometric flux and the new Hip- 
parcos parallax, we have determined the stellar radius and effective 
temperature nearly independent of theoretical atmospheric models. 
The constraints found in this work on the luminosity, T e « and ra- 
dius will be invaluable for critically testing theoretical pulsation 
models of a Cir. We will present this in our companion paper, 
where we report the detection of the large frequenc y spacing in 
a Cir , based on photometry from the WIRE satellite (j Bruntt et al.l 
l2008h . 

We analysed high-resolution spectra taken at two different ro- 
tational phases. The detailed abundance analyses yield the same 
mean abundances for the two spectra, but differences in the line 
shapes of certain elements give evidence for spots on the surface of 
a Cir. We confirm the a bundance pattern se en in a Cir in an earlier 
study using fewer lines (Kupk aetalJI 1996b . The results follow the 
general trend in other roAp stars, i.e. high photospheric abundance 
of Cr, Co, Y and the rare-earth elements Nd and Eu. 
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